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Abstract 

I An approach to quantum cosmology, relying on strengths of both canonical and path integral formalisms, is applied 

^ ' to the cosmological model, Bianchi type IX. Physical quantum states are constructed on the maximal slice of the 

' ' cosmological history. A path integral is derived which evolves observables off the maximal slice. This result is 

\ compared a path integral propagator derived earlier with conventional Faddeev-Poppov gauge fixing. 
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I. INTRODUCTION 
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Previously, we suggested a new approach to the quantum description of cosmology Relying on the strengths 
C ' of both canonical and path integral methods, this "composite approach" to quantization offers a real possibility 
, of calculating physical observables for cosmological models. The idea is to define the physical observables and inner 
product of the canonical theory and then use them to construct a path integral which evolves these physical observables. 
In the companion paper, we presented a path integral quantization of the cosmological model Bianchi type IX using 
O ■ Faddeev-Poppov techniques familiar from gauge theory. Here, we carry out this composite approach for the model 

cosmology Bianchi type IX and compare our result to that found previously . 
?-H One of the great difhculties in the canonical quantization of gravity is that, to find the physical quantities in a 

gauge invariant formalism, it is essentially necessary to solve the theory. However, as pointed out by Vince Moncrief 
^ , during a discussion three years ago at Santa Barbara, there are cases in which we can identify physical observables and 
define their algebra. For example, in spatially compact cosmological models, we can define the physical observables 
^ . to be the metric and extrinsic curvatures of the maximal slice, on which the trace of the extrinsic curvature vanishes. 
' As the metric and extrinsic curvature evaluated there label solutions uniquely, they coordinatize the physical phase 
space. Moreover, the Poisson brackets are the naive ones, evaluated on that surface. 

Alternatively, physical observables can be defined by a complete gauge fixing of the theory. At the classical level, 
the two procedures must agree. However, in the quantum theory, it is not obvious that either of these procedures 
must lead to a good quantum theory. Nor is it obvious that they will lead to the same quantum theory. Our main 
goal in this work is to examine and compare the quantum theories that follow from each of these procedures, for 
the Bianchi type IX model. In each case we construct the path integral formulation of the quantum theory from the 
associated canonical theory. We compare the resulting path integrals with each other as well as with the expressions 
we previously found using the Faddeev-Poppov procedure. 

We work with the Bianchi type IX cosmological model because it is the simplest model that is, as far as is known, 
not solvable in closed form. This means that, although the configuration space is two dimensional, it shares with 
the real theory the property that we know of no procedure to explicitly construct the solutions of the theory. The 
model describes a cosmological family with homogeneous but anisotropic spatial slices. These anisotropics are the 
two dynamical degrees of freedom. The model has been studied extensively, especially since Misner expressed the 
dynamics as single particle mechanics in a time dependent potential |^,^. A key property, which we shall exploit, is 
that all of this model's classical histories recoUapse [Q . Recently, quantizations of Bianchi type IX have been proposed 
by Kodama || (see also [|), Grahm §, and Marolf §. 

The paper begins with a quick review of the phase space of this model. Once the maximal slice is determined we 
present the physical quantization and path integral. The third section offers a comparison with our earlier work. The 
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paper finishes with discussion in Section 4. 

We present the derivation in "geometrized units," in which G = c = 1. 



II. SUMMARY OF THE CLASSICAL THEORY 



We consider class A Bianchi type IX models which describe anisotropic, homogeneous solutions to the vacuum 
Einstein equations.[j In Misner's chart on the phase space of Bianchi type IX the spatial metric is given by [|| 

h.,=e'-{e'^)^^ (1) 

where is diagonal and traceless matrix parameterized by diag(/3+ + \/3(3~ — \/3/3~, —2/3+) and a is the scale 
factor of the cosmology. As /3'^ is traceless, the scale factor for the cosmology, a, is related to the volume through 
Vh — e^". The action for the model Bianchi type IX in this chart is [|| 

/ = y p+dP+ + p^dp- + pada - \l^-Ne-^°"Hdt (2) 
in which H = \ (— +p+ + e''"Z^ .0 The triangularly shaped potential is given by 

U{f3^) = ie-«'3+ - cosh(2V3/3-) 

+ (cosh(4V3/3-) - l) . (3) 

Thus, mixmaster dynamics may be seen as the dynamics of a single particle in a time dependent potential. The action 
of Eq. (||) can be expressed by fixing the lapse, TV, so that N = .^y^^exp(3a). This allows us to write 

/ = / p+d/3+ + p-d(3^ + pada - HdX. (4) 



This action must be supplemented with the condition Ti, = 0. For a fixed value of this potential, the walls contract as 
the scale factor increases to the maximum volume slice. After this slice, the walls of the potential, for a fixed value, 
recede back from the center. 

To implement the proposed quantization of the theory, we need to determine the slice with vanishing extrinsic 
curvature. Since all Bianchi type IX classical histories end in recoUapse [Q, we can always identify this slice. The 
physical degrees of freedom are simply the two anisotropics, and In addition, on this slice the volume is 
maximal. Since -s/Zi = e^", to find this maximum volume slice, we must maximize a. From the action of Eq. (Q) we 
find, 

a = {a,H} ^ -Pa] (5) 

so that the volume is extremized (for finite volume) when pa vanishes. To find the maximum, note that 

aU=o = 37rN^e°'U{P^) (6) 

giving the maximum volume, e^", when the potential is less than zero. This defines an open region TZ, TZ :~ : 
Uif}^) < 0} shown in Fig. 1. This_ region is finite [p[. 

The kinematical phase space, F, is coordinatized by (a,/?^,/?^) and conjugate momenta {pa,p+,p~)- The 
physical phase space F is a four dimensional submanifold of F, which is defined as the quotient of F by the solution 
space of H = 0. We see that 



*The classification of Bianchi models concerns the irreducible parts of the structure constants of the isometry Lie group. 
Writing the structure constants in terms of these pieces, Cja' ~ ^jklS^^ + S^Vk], class A models are those for which Vi — 0. 

^Note that in the literature one often finds written V{I3±) := U(P±) + 1. This is convenient because V{I3±) is positive definite, 
however what is important to remember is that the actual potential U{f3±) is bounded from below by —1. 
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FIG. 1. The zero contour for the classical potential for Bianchi IX. Three "channels," one along the positive /? axis, and 
the other two sloping diagonally off in the negative /3~ direction extend to infinity. The area of the interior region, R, is finite. 



r = T*n (7) 

coordinatized by G TZ and p± . 

The coordinates of T correspond to the metric and extrinsic curvature of solutions on the unique maximal slice 
Pa — 0. The quantities O = {a,/9^,p±} on a slice, say pa = t, could be found by integrating the solutions from 
Pa — Q to Pa — T using the Hamilton's equation of motion (See, for example, We will denote this set of physical 

observables 0{t). As all solutions end in singularities in some finite time, the 0{t) are only defined on a subspace of 
r corresponding to data on the surface = which specifies non-singular solutions. 

For example, let us consider the observable for volume of the universe, on the slice where pa ~ t, found by solving 
the Hamiltonian constraint 



V{t) ^ e 



- p3a(r) _ 



(8) 



As each of the factors is a complicated function on F, the explicit form of this can can only be found by integration 
of the equation 

{F(p„),7^}-0 (9) 
on the kinematical phase space F with the initial condition, 

i.e. one must essentially solve the theory. At each r the subspace r(r) defined by V{t) > is the domain on which 
the C(t) are defined. That the physical observables are only defined on subspaces of f is a new feature of cosmology 
that raises questions such as whether is the Poisson bracket of these functions always defined. Of great interest is the 
extent to which this causes problems for the quantum theory, in either a canonical or a path integral formalism. 



III. MAXIMAL SLICE CANONICAL QUANTUM THEORY 

We first quantize the physical phase space F = T*TZ. Since the physical configuration variables, (/3^,/3^) £ 7?., 
and the conjugate momenta, (p+,p_) S , are coordinates on the classical phase space, take functions as 
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wavefunctions in the Hilbert space. The quantum theory is built from the configuration operators which act 
by multiplication, and the momentum operators pj. which act by differentiation pj. = —ihd± on the wavefunctions. 
These operators correspond to measurements of the corresponding quantities defined on the slice Pa — 0. Hence the 
subscript and superscript "0." 
We have the inner product, 

{X\0= J d^P^x(pf)af3^)- (11) 

K 

If the wavefunctions vanish on the boundary of TZ then and pj. are Hermitian in this inner product. 

Before considering the problem of observables at times other than pa = 0, it is interesting to study the operator 
that measures the volume on maximal slice. As this is a property of solutions, it is a physical observable, but as it is 
a function of data on the maximal slice it can be exhibited explicitly. It is convenient to consider the 4/3 power of 
the volume, which by Eq. (8) is 

written in terms of a metric defined by q^^ — q — l/{—\U\). This is non-singular in the interior but blows up 
on the boundary of TZ. Quantum mechanically, we want to represent V^^^ as an hermitian operator. This can be 
done by representing it as the Laplace-Beltrami operator for the metric qij on the region TZ. We thus defined the 
corresponding operator to be 

We note that the form of the inner product Eq. (11) may be interpreted to mean that the states are half-densities on 
TZ. 

As this is a hermitian operator acting on a finite area, we conjecture that the volume has a discrete spectra with 
finite degeneracy. I This is a familiar result of one dimensional quantum mechanics. 



IV. EVOLUTION IN GAUGE INVARIANT QUANTIZATION 



In principle, the entire physical quantum theory is contained in the states and operators on the physical Hilbert 
space. However, to make predictions we need to find observables (t) for nonzero r. If we were able to do so, all 
the physical information would be contained in A'^-time correlation functions of the form 

{Pf\d\r,)...d''{rj,)\(3f). (14) 

The problem is how to construct these operators given that we have neither their classical counterparts in closed form 
nor the proper subsets of F. As an example of such an operator, consider an operator V{t) corresponding to the 
volume of the universe on the slice where Pa = t. There are two approaches to this problem: 

(1.) In some appropriate approximation procedure, solve Eqs. (9) and (10) (or the appropriate equations for other 
observables) in the classical theory, and then, term by term in the approximation, define an ordering prescription which 
realizes the classical expression as a well-defined quantum mechanical operator. (We may note the approximation 
procedure must be able to keep track consistently of the regions of definition T{t) which satisfy T{ti) E r(r2) for 

T2 > Tl.) 

(2.) Evolve the operators quantum mechanically. As we have defined the quantum theory only for the physical 
observables, there is no operator that corresponds to Ti.; we don't have the quantum mechanical analogue of Eq. (9). 
We can associate a non-vanishing Hamiltonian to evolution in a physically meaningful time variable, and then use 
that to define an evolution operator. In the classical theory, this is straightforward. 

In the classical phase space, the Hamiltonian which evolves physical variables in pa is the conjugate quantity a, 
which is a time { Pa ~ t) dependent functional on F. This is given by 



"''This is delicate. On the boundaries of the non-compact region the metric factors, 1/ — \U\, diverge. The wavefunctions, 
interpreted as half-densities, must fall-off faster than the potential. 
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/i(t) = a(r) = 7 In 



P+(t)2+P_(t)2 



^21 



(15) 



On f , this satisfies {/i(t),p„} = 1> so that it generates evolution in on the constraint surface. We could construct 
an evolution operator in the physical Hilbert space if we could find the corresponding quantum operator. However, 
this is not easy as we do not know the operators $^{t) and p±(t); these themselves are supposed to be found by 
evolution. In some operator ordering, we would have to find operator solutions to the coupled operator equations. 



Mr) = -ln 



-\U0^)\ 



(16) 



and 



dT 



= [/3±(r),Mr)] 



(17) 



dp±{T) 
dT 



= [p±(T),/i(r)]. (18) 



Needless to say, as the theory cannot be solved at the classical level, this requires some approximation procedure. 
There arc also other issues associated with operator ordering ambiguities such as whether quantum observables should 
be hermitian, given that each classical universe has a finite lifetime. For, if we find a hermitian ordering for /i(t) then 
we ought to be able to evolve an arbitrary quantum state to arbitrarily large times. Given that these practical and 
conceptual difficulties face any attempt to proceed with this program, we turn to another approach to representing 
evolution based on the path integral. 



V. PATH INTEGRAL REPRESENTATION OF PHYSICAL EVOLUTION 

Path integration provides a practical resolution for these problems. While the observables 0{t) cannot be written 
in closed form on the phase space, they could be computed by summing over histories. That is, as a computer averages 
the effective action over an ensemble of paths generated by some statistical procedure, it can simply go to the points 
on the path where Pa = t and then tabulate the other observables at those events. In this way, averages such as Eq. 
(14) could be evaluated numerically, with an ensemble of paths with appropriate weights or measure. We study one 
such procedure in which the path integral measure may be derived from the physical quantum theory. In this way, 
the complementary strengths of the two approaches may be exploited, for a practical approach to calculating physical 
observables non-perturbatively in quantum cosmology. 

To construct the path integral corresponding to the matrix element 

{(if\0{T^)...0{T,Wf) (19) 

where r/ > r/ > Tj; / = 1, 2, 3, . . . , A^, we express the states in terms of physical states, elements of C^iTZ). These are 
the time evolved kets 

\PfTf)=U{TuTf)\Pfn) (20) 

where U (tj, t/) is the operator for evolution r. We derive a path integral expression for Eq.(19) in terms of a quantum 
effective action S = J^' L, a measure factor /x, and ranges of integration, 

{l3f\0{TN)...0{n)\pt) = [d'pd^pt,{p,p)] 0(T;v)...0(ri)e(^/'^)^ (21) 

where the 0{ti ) are the classical expressions, expressed as the value of each observable O on the path at the points 
Pa = T, and 

J\f= f [d^l3d^pnif3,p)]\e^'/''^^o (22) 
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(The brackets a notational device to indicate a factor of l/27r for each differential.) To do this we must first find a 
path integral expression for U{Tf,Ti). The hamiltonian that corresponds to the choice t = Pa [given in Eq. (16)], 
may be expressed as 



Mr) = 



In 



-|Zi(/3±)| 



In 



-p+ - pt 



+ 1 



The evolution operator is expressed in terms of the corresponding operator as 



U{0,t) = Texp 



^ Jo 



h{T')dT' 



(23) 



(24) 



Without the form of /3^(t) and P±{t), we do have the Hamiltonian /i(r) for finite r as an operator on the physical 
Hilbert space. To construct the path integral we make some conjectures about these operators and their spectra. 
These assumptions cannot be justified directly in the absence of further information about the solutions to Eqs. (16), 
(17), and (18), but they may be plausible in the light that the resulting path integral agrees with that constructed 
by the standard method of gauge fixed quantization. We assume that there exist simultaneous operator solutions to 
Eqs. (16), (17), and (18) such that they are all hermitian operators. As a result time evolution is unitary, so that the 
P'^ir) and p±{t) satisfy equal time commutation relations, [(3^ {t) , p±' (t)] = S±±'. In this case, for each r, there 
must be a complete basis of states made from eigenstates of ;5^(t). 
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/ d/3±|/3±(r))(/3±(r)| 



(25) 



where the range of the integral, 7?.(t), is the range of the spectra of the operators /3^(t). In the absence of a 
construction of the operators, we do not know this range, but we would like to argue that it is the whole of R"^. In 
the classical theory we have, from the Hamiltonian constraint, 



■4a 



(26) 



At r = 0, this restricts to the region TZ = TZ{/), However, there are initial classical configurations for all values of 
Pa] this means that for small r, the left hand side can be arbitrarily large, which means that there is no limit on how 
large the anisotropics can be at any finite t. As a result, if the theory is to have a good classical limit it seems 
likely that the spectra must be unbounded. 

Similarly, we assume that there is a complete set of states 

1= / d^f3\p±{T)){p±{T)\. (27) 

Given all these assumptions we may deduce 

(p±(T)|/3±(T))=e'('5^f++''"f-). (28) 

We use standard techniques to construct the path integral. Inserting complete sets of states in the usual convenient 
ways, we evaluate 



/i(r) = (/3±(r)|MT)b±(r)) 



1 



In 



-\U{P±)\ 



+ ln(-^±^ + l 



(29) 



to find that, 



f/(0,r) = j [d- 



^P'^d'^p± I exp 



u 

L 



P+I3++P-I3- -h{p±,(3^,T)dT 



(30) 



where the ranges of the integrals are, given our assumptions, unbounded. Comparing to Eq. (21), we then sec that in 
these coordinates the measure is trivial, the integration regions are unbounded and the effective lagrangian is given 
simply by L = p+/3+ + - h{p±,f3^,T). 
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VI. SCALE FACTOR TIME 



To compare this result to the earher derivation we'U change from an extrinsic curvature time gauge to a gauge 
related to the volume (r = InVTi = a). Geometrically, the measure transfers from one slice to another along the 
orbits of the gauge transformation. Unlike Pa, a is not monotonic on the whole history of the cosmology. When the 
gauge transformation is made we will choose which half of the history to consider. Expanding the phase space to 
include time and the Hamiltonian, so that 



;7(0,r) = J [(fl3'^(fp±dadpaS{T-pa)S{a-h{T))]exp 

The new gauge, r = a, may be introduced as 

1 = A„ y dA 5 (r - a^) 

where is the gauge transformed a. Under small gauge transformations, 

= a + X{a,n} + 0{X'^). 



- I p+P++p-f3 +PaadT 



(31) 



(32) 



(33) 



Thus, Aq = |{a, Ti.}\ = \pa\- Finally, we express the delta function involving a in terms of the Hamiltonian constraint 
as (5 (a — ^(t)) = SlH)\{H,Pa}\- Inserting these identities into the propagator of Eq. ( ^ gives. 



f7(0,T) = J [d^p^d^p±dadpc,dXS{T~pa)S{n)S{T~a^)\{a,n}\\{H,Pa}\] 



X exp 



-\p+P++P-l3 +paadT 



Observe that, as before with the new time choice, the extrinsic time choice satisfies 



dXS (r-p^) 



(34) 



(35) 



By explicitly changing gauge one can check that Ap^ is gauge independent and equals (for small gauge 

transformations). Performing an inverse gauge transformation to change gauge we have 



{7(0, r) 



d^ l3^d^p±dadpadX5 [t - pI ^\ 5 [U) 5 {t ^ a^) \{'H,a}\Ap, 



X exp 



- / P+I3+ +p^l3~ +PaadT 



(36) 



Performing the integration over the gauge parameter, allowing a delta to eat up the integration over a, and exponen- 
tiating the Hamiltonian constraint we find. 



U{0,t) = J [d^p^d^p±dpMdN] 

p+(i+ +P-I3- +pa~ NHdT 



X exp 



(37) 



the propagator in the physical phase space. This is the path integral for the canonical ADM Hamiltonian under 
the choice t = a. With this choice the Hamiltonian is given by —pa (the last term in Eq. (^)) and the lagrange 
multiplier or lapse is determined so that d = 1. The propagator of Eq. (|3^) can be integrated to a path integral in 
the phase space f ; as in j^] the integration over pa may be performed yielding 



UiO,T) 



d [i d P±^ 



exp 



- / P+P+ + p-P- - NHdT 



(38) 
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with F = [p^ + p?. + e4"W(/3±)] This can be directly compared with the propagator in Section III of ^ . These 
propagators are identical up to a measure factor fi{/3^). There is a simple reason for this: in the previous paper Q we 
began with a path integral which included an integration over all the components of the frame fields and connections, 
where as here we specified the theory in Eqs. (1-4) directly in terms of diagonal gauge. Thus, the path integral in 
[01 has an additional factor in the measure that came from the gauge fixing down to diagonal gauge. That factor is 
non-trivial, as can be seen directly from Eq. (40) of Had we begun the treatment of with the model defined 
in diagonal gauge, the results would have been identical to those found here. We conclude that the definition of the 
path integral depends on the point at which the reduced theory is taken to define the quantum theory; reduction does 
not commute with quantization. 

This measure factor is not relevant for the main question of interest here, which is whether this "composite" 
approach for defining the physical theory leads in the end to the same path integral as did the Faddeev-Poppov 
ansatz. The answer is that, at least in this case, if one starts from the same model, they are the same. 



VII. CONCLUSIONS 



We have implemented this new approach to quantum cosmology in the restricted, but non-trivial, case of the Bianchi 
type IX model. By first, identifying the physical degrees of freedom on a slice and then evolving physical observables 
off this slice by a sum over histories approach, we have shown how it is possible to learn physics without "solving the 
theory." We hope that this method will aid calculation of observables in other models and even the full theory. Given 
the physical theory on one slice, physical observables could be computed numerically. 

This work suggests several directions for further study. While the technology appears robust enough to handle the 
non-linear behavior of Bianchi type IX, it would be interesting if a full quantization could be carried out for one of 
the more simple Bianchi cosmological models. This composite approach also holds some hope for numerical work 
in that operators at arbitrary times may be found by computing the propagator for the cosmology. Expressed in 
terms of a path integral, this could be the basis for a numerical approach to quantum cosmology that has a physical 
interpretation grounded in the canonical theory. 
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